TITLE OF THE INVENTION 
APPARATUS AND METHOD FOR 
A RADIATION IMAGE THROUGH A GRID 

5 This is a continuation-in-part application of U.S. 

Patent Application No. 09/773, 095 filed on January 31, 2001 
entitled ''IMAGE ACQUISITION METHOD AND APPARATUS", 

FIELD OF THE INVENTION 
10 The present invention relates to an apparatus and a 

method for a radiation image of an object to be imaged 
through a grid for removing scattered rays from the object. 

BACKGROUND OF THE INVENTION 

15 With the recent trend toward the digitization of 

medical X-ray images. X-ray intensity spatial 
distributions can be acquired as digital images. For 
example^ currently available schemes include a scheme of 
forming a latent image on a photostimulable phosphor by 

20 using X-ray energy and acquiring an image from a laser 
pumping light distribution, a scheme of converting an X-ray 
intensity distribution into a light intensity distribution 
(fluorescence) , directly converting the distribution into 
an electrical signal by using a surface center having a 

25 plurality of pixels, and converting the signal into a 

digital image, and a scheme of directly converting an X-ray 
intensity distribution into a charge distribution. 
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The digitization of X-ray images has the following 
merits : 

• Storage and transfer can be efficiently performed. 

• Optimal images can be easily formed by digital image 
5 processing (a recovery from a failure in imaging operation 

can be made) . 

• Efficient diagnosis can be carried out. 

• A reduction in the cost of diagnosis can be attained. 
The problem of scattered X-rays produced when X-rays 

10 pass through an object has not been satisfactorily solved. 
To reduce the influences of scattered X-rays and obtain a 
high-contrast image, an optimal means is to use a scattered 
ray removing grid having many lead plates arranged in the 
same direction as in a conventional scheme using 

15 silver-halide films. 

Fig. 10 is a schematic sectional view of a structure 
using a grid- Reference numeral 81 denotes a point (X-ray 
focal point) of an X-ray tube from which X-rays are 
generated; 82, an object; 83, a grid; 84, an energy 

20 conversion unit for converting an X-ray intensity 

distribution into a light intensity or charge amount; and 
85, a sensor unit for spatially sampling the distribution. 
Scattered X-rays reflected by the object reach the grid 83 
as well as X-rays directly emitted from the X-ray tube. 

25 Most of the scattered X-rays are cut by the grid 83 made 
up of lead members facing the X-ray focal point 81. 

A drawback of the grid 83 is that it partly cuts direct 



- 2 - 



x-rays while it cuts scattered X-rays. This cutting 
pattern corresponds to the arrangement of the lead members 
of the grid, and the image generally suffers a stripe 
pattern. An X-ray image has evolved into (1) an image 
5 formed by a film-screen system (analog image) — ♦ (2) a 
digital image formed by reading a latent image formed by 
a photostimulable phosphor by laser scannings (3) an image 
formed by direct sampling (flat panel sensor) the 
two-dimensional spatial distribution of an X-ray dose in 
10 a two-dimensional space. Different measures have 

therefore been taken against a stripe pattern (grid image) 
on an image which originates from the lead members used for 
the grid 83. 

In the film-screen system in (1), the following two 
15 methods are available, which are used to remove a grid image 
or prevent interference with observation. 

(a) The grid itself is moved during radiation of 
X-rays to prevent the formation of a grid image while 
removing scattered rays. 
20 (b) The spatial frequency of the grid stripe pattern 

is increased to make it difficult for the human eye to 
perceive a grid image if it is formed on an image or prevent 
the grid image from overlapping the frequency component of 
image information . 
25 The means of moving the grid itself in (a) is effective 

in all cases of X-ray image acquisition. However, this 
means is difficult to use because of an increase in cost 
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due to, for example, a driving system for moving the grid, 
an increase in apparatus size, the relationship between the 
driving timing and the X-ray radiation timing, adjustment 
of the driving speed, and the like. 
5 The means of increasing the spatial frequency of the 

grid stripe pattern in (b) has its own limit. That is, when 
the frequency of the grid stripe pattern is set to a high 
spatial frequency at which no grid image is formed, since 
the thickness of each lead plate for blocking scattered rays 

10 is almost fixed, an area through which direct rays are 
transmitted narrows, and the use efficiency of the X-ray 
dose extremely decreases. As a consequence, imaging 
operation cannot be properly performed. 

In the era in which a latent image formed by a 

15 photostimulable phosphor was read by laser scanning and 
digitized, the idea of using anti-aliasing filter before 
sampling was introduced as a method of removing a grid image • 
When a latent image formed by a photostimulable phosphor 
is to be read by laser scanning and digitized, the image 

20 is scanned in a one-dimensional direction by a laser to 
temporarily form a signal form like a video signal, and the 
signal is sampled on the time axis. The frequency of the 
grid stripe pattern is increased to a certain degree, and 
laser scanning is performed in a direction perpendicular 

25 to the grid stripe pattern to form the grid stripe pattern 
into a periodic signal on the video signal. A grid image 
can be removed by the general idea of using an anti-aliasing 
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filter, i.e., performing sampling on the time axis after 
low-pass filtering in the state of an analog signal as this 
video signal. A similar method is disclosed in Japanese 
Patent No. 2507659, in which a grid image and its frequency 
5 are obtained by Fourier transformation of an image obtained 
by preliminary sampling, and a low-pass filter 
corresponding to the result is selected to remove a grid 
image . 

According to another method, sampling is performed 
10 on the time axis at intervals shorter than desired intervals 
instead of performing analog low-pass filtering to - 
eliminate aliasing of grid stripe pattern information, and 
the resultant information is separated from image 
information. Thereafter, digital low-pass filtering is 
15 performed, and the resultant image is digitally decimated 
(sub-sampled) , thereby obtaining an image at the desired 
sampling intervals. Similar methods are disclosed in 
Japanese Patent No. 2754068 and Japanese Patent Laid-Open 
No. 8-088765. 

20 In the advanced era in which a digital X-ray image 

can be obtained by directly sampling (using a flat panel 
sensor) a two-dimensional spatial distribution of an X-ray 
dose in (3) in a two-dimensional space, the above 
anti-aliasing filtering cannot be used. That is, a flat 

25 panel sensor is made up of a plurality of semiconductor 
pixels, and the two-dimensional spatial sampling pitch of 
the sensor cannot be reduced more than necessary in 



- 5 - 



consideration of technique and cost. The above idea of 
using an anti-aliasing filter cannot be applied to this 
method. The method disclosed in Japanese Patent Laid-Open 
No. 9-75332 is aimed at removing grid stripe pattern 
5 information in contrast to the method of obtaining an X-ray 
image by direct sampling in a two-dimensional space. In 
this method, the intervals of grid lead members are 
perfectly matched with the sampling pitch to match areas 
where direct X-rays are blocked by the grid stripe pattern 

10 with the gaps between the pixels, thereby preventing a grid 
stripe pattern from appearing on an image. 

Japanese Patent Laid-Open No. 9-98970 and 
U.S. Patent No. 5,801,385 disclose methods of setting the 
grid lead member intervals to be smaller than the sampling 

15 pitch and equal to or near the width of the opening of a 
light-receiving portion of one pixel, thereby reducing the 
contrast of a grid stripe pattern. In the method disclosed 
in U.S. Patent No. 5, 050, 198, a grid image is input and 
stored under a plurality of conditions. When imaging 

20 operation is actually performed by using the grid, the 
obtained image is divided by a grid image of the stored grid 
images which corresponds to the condition under which the 
actual imaging operation is performed, thereby removing the 
grid image. 

25 In the method disclosed in Japanese Patent Laid-Open 

No. 9-75332, which corresponds to the above technique of 
obtaining a digital X-ray image by direct sampling in a 
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two-dimensional space using a flat panel in the 
two-dimensional space, it is very difficult to perfectly 
match the grid lead member intervals with the sampling pitch. 
In the methods disclosed in Japanese Patent Laid-Open 
5 No. 9-98970 and U. S. PatentNo. 5, 801, 385, a grid image can 
be effectively removed by reducing the grid lead member 
intervals below the sampling pitch to be equal or near the 
width of the opening of a light-receiving portion of one 
pixel. However, as the flat panel sensor increases in 

10 resolution, and the sampling pitch becomes 0.1 mm or less, 
the grid lead member intervals are required to be very 
small; 10 or more grid lead members per mm. If the intervals 
become so small, since the thickness of each lead plate for 
blocking scattered rays is almost fixed, the areas through 

15 which direct rays pass narrow, and the use efficiency of 
the X-ray dose becomes extremely low. As a consequence, 
proper imaging operation cannot be performed. 

SUMMARY OF THE INVENTION 
20 It is an object of the present invention to acquire 

a radiation image in which the interference with 
observation due to a stripe pattern , which originates from 
a scattered ray removing grid, is suppressed. 

It is another object of the present invention to remove 
25 a stripe pattern, which originates from a scattered ray 
removing grid, in a radiation image, as needed. 

According to an aspect of the present invention, there 



- 7 - 



• 'I 



is provided a radiation image acquisition apparatus 
characterized by comprising a sensor for spatially sampling 
a radiation transmission distribution of an object to be 
imaged at a spatial sampling interval and acquiring an image 
5 of said object, and a scattered ray removing grid for 
removing scattered rays from said object, wherein an 
interval of elements of said scattered ray removing grid 
is set such that a spatial frequency of a stripe pattern, 
in said image, which originates from said scattered ray 

10 removing grid becomes not less than 1/3 and not greater than 
4 0% of a sampling frequency that is a reciprocal of said 
spatial sampling interval. 

It is another aspect of the present invention, there 
is provided a radiation image acquisition method of 

15 spatially sampling a radiation transmission distribution 
of an object to be imaged at a spatial sampling interval 
and acquiring an image of said object using a sensor and 
a scattered ray removing grid for removing scattered rays 
from said object characterized by comprising setting an 

20 interval of elements of said scattered ray removing grid 
such that a spatial frequency of a stripe pattern, in said 
image, which originates from said scattered ray removing 
grid becomes not less than 1/3 and not greater than 40% of 
a sampling frequency that is a reciprocal of said spatial 

25 sampling interval. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 
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apparatus characterized by comprising a sensor for 
spatially sampling a radiation transmission distribution 
of an object to be imaged at a spatial sampling interval 
and acquiring an image of said object, and a scattered ray 
5 removing grid for removing scattered rays from said object, 
wherein an interval of elements of said scattered ray 
removing grid is set such that a spatial frequency that is 
a reciprocal of said interval of elements of said scattered 
ray removing grid is within fs(n+l/3) fs(n+0.4) or 

10 fs{n+0.6) fs(n+2/3), where 1/fs is said spatial sampling 
interval and n is an integer not less than 0. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 
method of spatially sampling a radiation transmission 

15 distribution of an object to be imaged at a spatial sampling 
interval and acquiring an image of said object using a 
sensor and a scattered ray removing grid for removing 
scattered rays from said object characterized by comprising 
setting an interval of elements of said scattered ray 

20 removing grid such that a spatial frequency that is a 

reciprocal of said interval of elements of said scattered 
ray removing grid is within fs(n+l/3) fs(n+0.4) or 
fs(n+0.6) fs(n+2/3), where 1/fs is said spatial sampling 
interval and n is an integer not less than 0. 

25 According to still another aspect of the present 

invention, there is provided a radiation image acquisition 
apparatus characterized by comprising an image acquisition 
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unit for spatially sampling a radiation transmission 
distribution of an object to be imaged through a scattered 
ray removing grid for removing scattered rays from said 
object at a spatial sampling interval and acquiring an image 
5 of said object, an image processing unit for removing a 
stripe pattern originating from said scattered ray removing 
grid from said image by image processing, and a selection 
unit for allowing selection between removal and nonremoval 
of said stripe pattern by using said image processing unit, 

10 wherein the selection is performed in accordance with an 
application purpose. 

According to still another aspect of the present 
invention^ there is provided a radiation image acquisition 
apparatus characterized by comprising an image acquisition 

15 unit for spatially sampling a radiation transmission 

distribution of an object to be imaged through a scattered 
ray removing grid for removing scattered rays from said 
object at a spatial sampling interval and acquiring an image 
of said object, an image processing unit for removing a 

20 stripe pattern originating from said scattered ray removing 
grid from said image by image processing, and a selection 
unit for allowing selection between removal and nonremoval 
of said stripe pattern by using said image processing unit, 
wherein the selection is performed in accordance with a 

25 portion to be imaged of said object. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 
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apparatus characterized by comprising an image acquisition 
unit for spatially sampling a radiation transmission 
distribution of an object to be imaged through a scattered 
ray removing grid for removing scattered rays from said 
5 object at a spatial sampling interval and acquiring an image 
of said object, an image processing unit for removing a 
stripe pattern originating from said scattered ray removing 
grid from said image by image processing, and a selection 
unit for allowing selection between removal and nonremoval 

10 of said stripe pattern by using said image processing unit, 
wherein the selection is performed in accordance with an 
amplitude of said stripe pattern. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 

15 apparatus characterized by comprising an image acquisition 
unit for spatially sampling a radiation transmission 
distribution on an object to be imaged through a scattered 
ray removing grid for removing scattered rays from said 
object at a spatial sampling interval and acquiring an image 

20 of said object, an image processing unit for removing a 
stripe pattern originating from said scattered ray removing 
grid from said image by image processing, and a selection 
unit for allowing selection between removal and nonremoval 
of said stripe pattern by using said image processing unit, 

25 wherein the selection is performed in accordance with a 
magnitude of contrast of said stripe pattern. 

According to still another aspect of the present 
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invention, there is provided a radiation image acquisition 
method characterized by comprising a step of spatially 
sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 
5 removing scattered rays from said object at a spatial 
sampling interval and acquiring an image of said object, 
a step of selecting between removal and nonremoval of a 
stripe pattern originating from said scattered ray removing 
grid from said image by image processing, wherein the 

10 selection is performed in accordance with an application 
purpose, and a step of removing said stripe pattern by image 
processing in accordance with the selection. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 

15 method characterized by comprising a step of spatially 
sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
sampling interval and acquiring an image of said object, 

20 a step of selecting between removal and nonremoval of a 
stripe pattern originating from said scattered ray removing 
grid from said image by image processing, wherein the 
selection is performed in accordance with a portion to be 
imaged of said object, and a step of removing said stripe 

25 pattern by image processing in accordance with the 
selection . 

According to still another aspect of the present 
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invention, there is provided a radiation image acquisition 
method characterized by comprising a step of spatially 
sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 
5 removing scattered rays from said object at a spatial 
sampling interval and acquiring an image of said object, 
a step of selecting between removal and nonremoval of a 
stripe pattern originating from said scattered ray removing 
grid from said image by image processing, wherein the 
10 selection is performed in accordance with an amplitude of 
said stripe pattern, and a step of removing said stripe 
pattern by image processing in accordance with the 
selection. 

According to still another aspect of the present 
15 invention, there is provided a radiation image acquisition 
method characterized by comprising a step of spatially 
sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
20 sampling interval and acquiring an image of said object, 
a step of selecting between removal and nonremoval of a 
stripe pattern originating from said scattered ray removing 
grid from said image by image processing, wherein the 
selection is performed in accordance with a magnitude^ of 
25 contrast of said stripe pattern, and a step of removing said 
stripe pattern by image processing in accordance with the 
selection. 
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According to still another aspect of the present 
invention^ there is provided a radiation image processing 
apparatus for processing an image acquired by spatially 
sampling a radiation transmission distribution of an object 
5 to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
sampling interval characterized by comprising an image 
processing unit for removing a stripe pattern originating 
from said scattered ray removing grid from said image by 
10 image processing, and a selection unit for allowing 

selection between removal and nonremoval of said stripe 
pattern by using said image processing unit, wherein the 
selection is performed in accordance with an application 
purpose. 

15 According to still another aspect of the present 

invention, there is provided a radiation image processing 
apparatus for processing an image acquired by spatially 
sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 

20 removing scattered rays from said object at a spatial 
sampling interval characterized by comprising an image 
processing unit for removing a stripe pattern originating 
from said scattered ray removing grid from said image by 
image processing, and a selection unit for allowing 

25 selection between removal and nonremoval of said stripe 
pattern by using said image processing unit, wherein the 
selection is performed in accordance with a portion to be 
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imaged of said object. 

According to still another aspect of the present 
invention, there is provided a radiation image processing 
apparatus for processing an image acquired by spatially 
5 sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
sampling interval characterized by comprising an image 
processing unit for removing a stripe pattern originating 

10 from said scattered ray removing grid from said image by 
image processing, and a selection unit for allowing- 
selection between removal and nonremoval of said stripe 
pattern by using said image processing unit, wherein the 
selection is performed in accordance with an amplitude of 

15 said stripe pattern. 

According to still another aspect of the present 
invention, there is provided a radiation image processing 
apparatus for processing an image acquired by spatially 
sampling a radiation transmission distribution of an object 

20 to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
sampling interval characterized by comprising an image 
processing unit for removing a stripe pattern originating 
from said scattered ray removing grid from said image by 

25 image processing, and a selection unit for allowing 

selection between removal and nonremoval of said stripe 
pattern by using said image processing unit, wherein the 
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selection is performed in accordance with a magnitude of 
contrast of said stripe pattern. 

According to still another aspect of the present 
invention, there is provided a radiation image processing 
5 method of processing an image acquired by spatially 

sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
sampling interval characterized by comprising a step of 

10 selecting between removal and nonremoval of a stripe 

pattern originating from said scattered ray removing grid 
from said image by image processing, wherein the selection 
is performed in accordance with an application purpose, and 
a step of removing said stripe pattern by image processing 

15 in accordance with the selection. 

According to still another aspect of the present 
invention, there is provided a radiation image processing 
method of processing an image acquired by spatially 
sampling a radiation transmission distribution of an object 

20 to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
sampling interval characterized by comprising a step of 
selecting between removal and nonremoval of a stripe 
pattern originating from said scattered ray removing grid 

25 from said image by image processing, wherein the selection 
is performed in accordance with a portion to be imaged of 
said object, and a step of removing said stripe pattern by 
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image processing in accordance with the selection. 

According to still another aspect of the present 
invention^ there is provided a radiation image processing 
method of processing an image acquired by spatially 
5 sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 
removing scattered rays from said object at a spatial 
sampling interval characterized by comprising a step of 
selecting between removal and nonremoval of a stripe 

10 pattern originating from said scattered ray removing grid 
from said image by image processing^ wherein the selection 
is performed in accordance with an amplitude of said stripe 
pattern, and a step of removing said stripe pattern by image 
processing in accordance with the selection. 

15 According to still another aspect of the present 

invention, there is provided a radiation image processing 
method of processing an image acquired by spatially 
sampling a radiation transmission distribution of an object 
to be imaged through a scattered ray removing grid for 

20. removing scattered rays from said object at a spatial 

sampling interval characterized by comprising a step of 
selecting between removal and nonremoval of a stripe 
pattern originating from said scattered ray removing grid 
from said image by image processing, wherein the selection 

25 is performed in accordance with a magnitude of contrast of 
said stripe pattern, and a step of removing said stripe 
pattern by image processing in accordance with the 
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selection. 

According to still another aspect of the present 
invention, there is provided a design method of designing 
at least one of a sensor and a scattered ray removing grid 
5 used for an apparatus for spatially sampling a radiation 
transmission distribution of an object to be imaged at a 
spatial sampling interval and acquiring an image of said 
object using said sensor and said scattered ray removing 
grid for removing scattered rays from said object 

10 characterized by comprising determining at least one of an 
interval of elements of said scattered ray removing grid 
and said spatial sampling interval of said sensor such that 
a spatial frequency of a stripe pattern, in said image, 
which originates from said scattered ray removing grid 

15 becomes not less than 1/3 and not greater than 40% of a 
sampling frequency that is a reciprocal of said spatial 
sampling interval. 

According to still another aspect of the present 
invention, there is provided a design method of designing 

20 at least one of a sensor and a scattered ray removing grid 
used for an apparatus for spatially sampling a radiation 
transmission distribution of an object to be imaged at a 
spatial sampling interval and acquiring an image of said 
object using said sensor and said scattered ray removing 

25 grid for removing scattered rays from said object 

characterized by comprising determining at least one of an 
interval of elements of said scattered ray removing grid 
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and said spatial sampling interval of said sensor such that 
a spatial frequency that is a reciprocal of said interval 
of elements of said scattered ray removing grid is within 
fs(n+l/3) - fs(n+0.4) or fs(n+0.6) - fs(n+2/3), where 1/fs 
5 is said spatial sampling interval and n is an integer not 
less than 0. 

According to still another aspect of the present 
invention, there is provided a radiographic apparatus 
characterized by comprising a sensor for spatially sampling 

10 a radiation transmission distribution of an object to be 
imaged at a spatial sampling interval and acquiring an image 
of said object/ and a grid for reducing scattered radiation 
from said object , wherein a sampling frequency Fs of said 
sensor as a reciprocal of the spatial sampling interval and 

15 a spatial frequency Fg of said grid as a reciprocal of an 
interval of shades of elements of said grid on an 
image-receiving surface of said sensor substantially 
satisfy Fg = j*Fs/3, where j is a positive integer except 
for multiples of three. 

20 According to still another aspect of the present 

invention, there is provided a radiographic method of 
spatially sampling a radiation transmission distribution 
of an object to be imaged at a spatial sampling interval 
by using a sensor and a grid for reducing scattered 

25 radiation from said object, and acquiring an image of said 
object characterized by comprising acquiring an image of 
said object such that a sampling frequency Fs of said sensor 

- 19 - 



as a reciprocal of the spatial sampling interval and a 
spatial frequency Fg of said grid as a reciprocal of an 
interval of shades of elements of said grid on an 
image-receiving surface of said sensor substantially 
5 satisfy Fg = j-Fs/3, where j is a positive integer except 
for multiples of three. 

According to still another aspect of the present 
invention, there is provided a design method of designing 
at least one of a sensor and a grid used for a radiographic 

10 apparatus for spatially sampling a radiation transmission 
distribution of an object to be imaged at a spatial sampling 
interval and acquiring an image of said object using said 
sensor and said grid for reducing scattered rays from said 
object characterized by comprising determining at least one 

15 of the spatial sampling interval of said sensor and an 
interval of elements of said grid such that a sampling 
frequency Fs of said sensor as a reciprocal of the spatial 
sampling interval and a spatial frequency Fg of said grid 
as a reciprocal of an interval of shades of elements of said 

20 grid on an image-receiving surface of said sensor 

substantially satisfy Fg = j*Fs/3, where j is a positive 
integer except for multiples of three. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 

25 apparatus characterized by comprising a sensor for 

spatially sampling a radiation transmission distribution 
of an object to be imaged at a spatial sampling interval 
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and acquiring an image of said object, and a scattered ray 
removing grid for removing scattered rays from said object, 
wherein an interval of elements of said scattered ray 
removing grid is set such that a spatial frequency of a 
5 stripe pattern, in said image, which originates from said 
scattered ray removing grid becomes not greater than 40% 
of a sampling frequency that is a reciprocal of said spatial 
sampling interval . 

According to still another aspect of the present 

10 invention, there is provided a radiation image acquisition 
method of spatially sampling a radiation transmission 
distribution of an object to be imaged at a spatial sampling 
interval and acquiring an image of said object using a 
sensor and a scattered ray removing grid for removing 

15 scattered rays from said object characterized by comprising 
setting an interval of elements of said scattered ray 
removing grid such that a spatial frequency of a stripe 
pattern, in said image, which originates from said 
scattered ray removing grid becomes not greater than 40% 

20 of a sampling frequency that is a reciprocal of said spatial 
sampling interval . 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 
apparatus characterized by comprising a sensor for 

25 spatially sampling a radiation transmission distribution 
of an object to be imaged at a spatial sampling interval 
and acquiring an image of said object, and a scattered ray 
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removing grid for removing scattered rays from said object, 
wherein an interval of elements of said scattered ray 
removing grid is set such that a spatial frequency that is 
a reciprocal of said interval of elements of said scattered 
5 ray removing grid is within fs(n+0.25) - fs(n+0.4) or 
f s (n+O . 6) f s (n+0 . 75) , where 1/f s is said spatial sampling 
interval and n is an integer not less than 0. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 

10 apparatus characterized by comprising a sensor for 

spatially sampling a radiation transmission distribution 
of an object to be imaged at a spatial sampling interval 
and acquiring an image of said object, and a scattered ray 
removing grid for removing scattered rays from said object, 

15 wherein an interval of elements of said scattered ray 

removing grid is set such that a spatial frequency that is 
a reciprocal of said interval of elements of said scattered 
ray removing grid is within fs(n+0.3) fs(n+0.4) or 
fs(n+0-6) - fs(n+0.7), where 1/fs is said spatial sampling 

20 interval and n is an integer not less than 0. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 
method of spatially sampling a radiation transmission 
distribution of an object to be imaged at a spatial sampling 

25 interval and acquiring an image of said object using a 
sensor and a scattered ray removing grid for removing 
scattered rays from said object characterized by comprising 
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setting an interval of elements of said scattered ray 
removing grid such that a spatial frequency that is a 
reciprocal of said interval of elements of said scattered 
ray removing grid is within fs(n+0.25) - fs(n+0.4) or 
5 fs(n+0.6) f s (n+0 . 75) , where 1/fs is said spatial sampling 
interval and n is an integer not less than 0. 

According to still another aspect of the present 
invention, there is provided a radiation image acquisition 
method of spatially sampling a radiation transmission 

10 distribution of an object to be imaged at a spatial sampling 
interval and acquiring an image of said object using a 
sensor and a scattered ray removing grid for removing 
scattered rays from said object characterized by comprising 
setting an interval of elements of said scattered ray 

15 removing grid such that a spatial frequency that is a 

reciprocal of said interval of elements of said scattered 
ray removing grid is within fs(n+0.3) fs(n+0.4) or 
f s (n+0 .6) f s (n+0 . 7) , where 1/fs is said spatial sampling 
interval and n is an integer not less than 0. 

20 According to still another aspect of the present 

invention, there is provided a design method of designing 
at least one of a sensor and a scattered ray removing grid 
used for an apparatus for spatially sampling a radiation 
transmission distribution of an object to be imaged at a 

25 spatial sampling interval and acquiring an image of said 
object using said sensor and said scattered ray removing 
grid for removing scattered rays from said object 
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characterized by comprising determining at least one of an 
interval of elements of said scattered ray removing grid 
and said spatial sampling interval of said sensor such that 
a spatial frequency of a stripe pattern, in said image, 
5 which originates from said scattered ray removing grid 
becomes not greater than 40% of a sampling frequency that 
is a reciprocal of said spatial sampling interval. 

According to still another aspect of the present 
invention, there is provided a design method of designing 

10 at least one of a sensor and a scattered ray removing grid 
used for an apparatus for spatially sampling a radiation 
transmission distribution of an object to be imaged at a 
spatial sampling interval and acquiring an image of said 
object using said sensor and said scattered ray removing 

15 grid for removing scattered rays from said object 

characterized by comprising determining at least one of an 
interval of elements of said scattered ray removing grid 
and said spatial sampling interval of said sensor such that 
a spatial frequency that is a reciprocal of said interval 

20 of elements of said scattered ray removing grid is within 
fs(n+0.25) - fs(n+0.4) or fs(n+0.6) - f s (n+0 • 75) , where 
1/fs is said spatial sampling interval and n is an integer 
not less than 0. 

According to still another aspect of the present 

25 invention, there is provided a design method of designing 
at least one of a sensor and a scattered ray removing grid 
used for an apparatus for spatially sampling a radiation 
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transmission distribution of an object to be imaged at a 
spatial sampling interval and acquiring an image of said 
object using said sensor and said scattered ray removing 
grid for removing scattered rays from said object 
5 characterized by comprising determining at least one of an 
interval of elements of said scattered ray removing grid 
and said spatial sampling interval of said sensor such that 
a spatial frequency that is a reciprocal of said interval 
of elements of said scattered ray removing grid is within 
10 fs(n+0.3) - fs(n+0.4) or fs(n+0.6) fs(n4-0.7) , where 1/fs 
is said spatial sampling interval and n is an integer not 
less than 0. 

Other objects and advantages besides those discussed 
15 above shall be apparent to those skilled in the art from 
the description of a preferred embodiment of the invention 
which follows. In the description, reference is made to 
accompanying drawings, which form a part thereof, and which 
illustrate an example of the invention. Such example, 
20 however, is not exhaustive of the various embodiments of 
the invention, and therefore reference is made to the claims 
which follow the description for determining the scope of 
the invention. 

25 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram for explaining the first 
embodiment of the present invention; 
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Figs. 2A and 2B are graphs for explaining spatial 
spectra in an image; 

Fig. 3 is a flow chart for the implementation of the 
first embodiment by software; 
5 Fig. 4 is a block diagram for explaining the second 

embodiment of the present invention; 

Fig. 5 is a block diagram for explaining the third 
embodiment of the present invention; 

Fig. 6 is a block diagram for explaining the fourth 
10 embodiment of the present invention; 

Figs. 7{l-a) to 7(3-b) are graphs for explaining 
artifacts originating from a grid stripe pattern; 

Figs. 8A and 8B are graphs for explaining beat noise 
produced when a grid stripe pattern is sampled; 
15 Fig. 9 is a graph showing an example of the spatial 

frequency range of a grid body; 

Fig. 10 is a view for explaining a scattered ray 
removing grid; 

Fig. 11 is a block diagram for explaining the fifth 
20 embodiment of the present invention; 

Fig. 12 is a view showing an example of a flat panel 
sensor; 

Fig. 13 is a graph showing the MTF of the flat panel 
sensor, the spectrum of a grid stripe pattern, and the like; 
25 Fig. 14 is a view showing the positional relationship 

between an X-ray source, grid, and flat panel sensor; and 

Fig. 15 is a graph for explaining sampling of the 
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spectra of the harmonic components of a grid stripe pattern. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The embodiments of the present invention will be 
5 described below with reference to the accompanying 
drawings . 

(First Embodiment) 

In an image acquisition apparatus according to the 
first embodiment of the present invention, when an acquired 
10 image is to be reproduced, the problem in the prior art is 
solved by setting the intervals of grid lead members such 
that even if a grid image exists and interferes with an image 
component to a certain degree^ the produced stripes have 
a frequency which hardly makes an observer have a sense of 
15 incongruity. 

This embodiment will be described below. In digital 
images obtained by sampling, image artifacts produced by 
the grid are classified into the following two categories: 
(A) grid images having low spatial frequencies; and 
20 (B) beat images obtained by sampling grid images 

having low spatial frequencies. 

Figs. 7(l-a) to7(3-b) schematically show the states 
of artifacts to explain artifacts belonging to categories 
(A) and (B) . For the sake of descriptive convenience, 
25 consider only a fundamental wavelength, assuming that the 
second- and higher-order harmonics of a grid stripe pattern 
are not resolved. 
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Fig. 7 (1-a) shows the state of an artifact belonging 
to category (A) in the frequency domain (only the positive 
region) . Reference numeral 61 denotes an image component 
region, which exhibits a substantially maximum frequency 
5 fi; and 62, an image harmonic component mathematically 
produced by sampling. In this case, a sampling frequency 
is represented by fs, and Nyquist frequency fnq = fs/2. A 
line spectrum 63 is a grid stripe pattern component having 
a frequency fg. As shown in Fig. 7 (1-a), since this 
10 component overlaps the image component 61 to impair the 
image quality, and has a low frequency, the observer- 
experiences a sense of incongruity from the stripe pattern. 

Fig. 7 (1-b) one-dimensionally shows how sampling is 
performed. In Fig. (1-b), (bullet) indicates a 

15 sampling point. 

It is obvious from the above description that the 
frequency fg of the grid stripe pattern and the frequency 
fi (fi < fs/2) of the image must satisfy: 

fg > i ... (1) 

20 Fig. 7(2-a) shows the state of another artifact 

belonging to category (A) in the frequency domain (only in 
the positive region) , A region 61 is an image component 
region and exhibits the maximum frequency fi. Since the 
frequency fg of the grid is higher than the Nyquist 
25 frequency fnq that is 1/2 the sampling frequency fs, a 
frequency component 64 of fs - fg, which is aliasing, 
appears. In the case shown in Fig. 7 (2-a) , the frequency 
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component 64 as aliasing overlaps the image component 61 
and impairs the image quality, and hence a stripe pattern 
gives the observer a sense of incongruity. Fig, 7 (2-b) 
one-dimensionally shows how sampling is performed. In 
5 Fig. (2-b), (bullet) indicates a sampling point, and 

the dashed line indicates a signal form before sampling. 
That is, to prevent impairment of the image component, the 
frequency fg of the grid stripe pattern and the frequency 
fi (fi < fs/2) of the image must satisfy: 

10 fs - fg > fi ... (2) 

According to inequalities (1) and (2), the frequency 
fg of the grid stripe pattern must satisfy the following 
inequality, in relation to the maximum frequency f i of the 
image and the sampling frequency fs: 

15 fi < fg < fs - fi (fi < fs/2) -..(3) 

In this case as well, an artifact belonging to 
category (B) may be produced. 

Fig. 7 (3-a) shows a state of an artifact when the grid 
stripe pattern frequency is set in the spatial frequency 

20 region which satisfies inequality (3). Fig. 7(3-b) 

one-dimensionally shows how sampling is performed. In 
Fig. (3-b) , "®" (bullet) indicates a sampling point, and 
the dashed line indicates a signal form before sampling. 
Referring to Fig. 7 (3-a), it seems that since an 

25 image frequency component 61 does not overlap a grid image 
component 63, the image is not impaired, and the observer 
has no sense of incongruity. If, however, beat noise-like 
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amplitude variation component overlaps a frequency 
component of image information as shown in Fig. 7 (3-b) , the 
observer of the image recognizes this overlap as an artifact, 
and a stripe pattern gives the observer a sense of 
5 incongruity. The frequency of this amplitude variation is 
given by | fs/2 - f g | . When some nonlinear conversion is 
performed for the image, this variation component may 
actually become an artifact that impairs the image. If, 
for example, the frequency of a grid stripe pattern on an 

10 image is set in a frequency region where beat-like artifacts 
(moire) occur (i.e., a high-frequency region), the grid 
stripe pattern and beat-like artifacts may be removed by 
filtering without damaging the image information (e.g., 
filtering to remove only components near the spatial 

15 frequency of the grid stripe pattern) . If, however, 
nonlinear element conversion (e.g., logarithmic 
conversion) is performed for a low-frequency beat-like 
artifact, the artifact becomes a signal component having 
an actual spectrum. In this case, this artifact cannot be 

20 removed by the above filtering technique. A beat-like 
artifact is repetitions of high-amplitude and low- (or 
zero-) amplitude portions of a stripe pattern. Since a 
high-amplitude portion inevitably differs in noise 
characteristics from a low-amplitude portion, even if 

25 stripe components are removed, variations in noise 

characteristics remain. As described above, even if the 
frequency of a grid stripe pattern on an image is set in 
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a frequency region where beat-like artifacts occur, and the 
grid stripe pattern is removed by filtering, some kind of 
mark of a beat-like artifact may remain. Therefore, the 
beat-like artifact itself (the occurrence or magnitude of 
5 the beat-like artifact) should be suppressed. 

Ideally, all artifact frequencies should fall 
outside the frequency fi of the image. In this case, no 
problem arises . Although the frequency f i is the maximum 
frequency of the image, this frequency may be regarded as 

10 a frequency that gives the observer no sense of incongruity 
when the image is reproduced (hard copy, monitor display, 
or the like) . Both inequality (3) and | fs/2 - f g | > fi 
must be satisfied at once. Assuming that fg is set at a 
position corresponding to fs > fg > fs/2, the following 

15 inequality is established: 

fi + fs/2 < fg < fs - fi ...(4) 
To satisfy this inequality, f i < f s/4 must also be satisfied, 
which is a strict condition. If fi < f s/4 is not satisfied 
in inequality (4), there is no overlap between the ranges 

20 defined by the left- and right-side signs. That is, an 
ideal condition for eliminating the influences of the grid 
is that the frequency band of the image is equal to or less 
than 1/2 the Nyquist frequency, and the frequency of the 
grid stripe pattern is near 1/2 the Nyquist frequency. This 

25 condition means that if an image to be acquired is 

determined, the sampling frequency must be set to at least 
four times higher than the maximum spatial frequency of the 
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image to be acquired. 

No consideration, however, is given to the power of 
beat noise-like variation component. The present inventor 
has contrived and proven that the condition given by 
5 inequality (4) can be moderated, on the basis of the result 
obtained by comparing the power of the beat noise-like 
variation component with the power of the grid stripe 
pattern and the result obtained by observing the image 
actually formed by using the grid. 

10 As the frequency of the stripe pattern separates 

farther from the Nyquist frequency, the amplitude variation 
as beat noise increases in frequency. As a consequence, 
this component cannot be observed. Consider how much this 
stripe pattern frequency differs from the Nyquist frequency 

15 when it becomes difficult for the observer to observe the 
variation component. 

Line spectra always exist at positions which have a 
mirror-image relationship centering on the Nyquist 
frequency whenever sampling is performed, and beat noise 

20 is always produced between them. In the cases shown in 
Figs. 7(l-b) and 7(2-b), although beat noise is produced, 
it does not exist apparently. One conceivable reason for 
this is attributed to frequency; in the case shown in 
Fig. 7 (2-a) or 7 (3-a) , the distance (frequency) between the 

25 two spectra having a mirror-image relationship centering 
on the Nyquist frequency is large, and hence the beat noise 
frequency becomes high. At this time, the spectrum of the 
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fundamental sine wave exists sufficiently below the Nyquist 
frequency, and they greatly differ in frequency and power. 
For this reason, only the fundamental sine wave having 
stronger power is strongly recognized by the observer. 
5 The manner in which beat noise is produced is 

mathematically expressed. Consider a case wherein a sine 
wave having the spatial frequency fg is sampled at the 
sampling frequency fs. Assume that fg > fs/2 is set in 
consideration of the grid actually used. This is not a 

10 necessary condition. In this case, as shown in Fig. 8A, 
a line spectrum pair corresponding to the frequency of the 
sine wave is produced. Let a/2 be the peak of each line 
spectrum, and a x cos(2 7cfgx) be the initial grid image. 
When a cosine wave having the frequency fg is sampled at 

15 the frequency fs, two cosine waves appear below fs. 

g(x) = a{cos (2 7c f gx) + cos(27c(fs - fg)x)} = 2a x 
cos (27r(fg - fs/2)x) x cos(27csx/2) ...(5) 
Equation (5) represents beat noise, which is 
equivalent to the amplitude modulated by a sine wave having 

20 a frequency corresponding to the difference between the two 
sine waves. As (fg/2 - fs) becomes a small value other than 
0, an unstable amplitude variation (beat) with a low 
frequency occurs. 

According to the Shannon's sampling theorem, data 

25 sampled at a frequency equal to or lower than the Nyquist 
frequency can be completely reconstructed by using an ideal 
filtering means (a filter that passes signals having 



- 33 - 



frequencies equal to or lower than the Nyquist frequency) , 
and no beat is produced by the line spectrum pair. It is 
likely that the observer strongly recognizes the beat noise 
in the case shown in Fig, 7(3-b) because of the filtering 
5 means . 

In general, a sampled signal is reconstructed by 
connecting the sampling points with a straight line or the 
like. This means differs from an ideal filtering means 
according to the sampling theorem. 
10 A person with normal visual perception or display 

apparatus interpolates by connecting neighboring points 
with a straight line without using any ideal filter 
(convolution using a sine function as a kernel) as in the 
sampling theorem. The same applies to a case wherein a 
15 signal is observed as an image. That is, such a difference 
between normal visual perception and the sampling theorem 
appears as beat noise. 

Interpolating with a straight line amounts to 
filtering with a characteristic like a characteristic curve 
20 71 in Fig. 8B. 

According to the characteristic curve 71 in Fig. 8B, 
a form s(f) of the filter is given by 

s(f) = sin^(7cf/fs)/(7rf/fs)2 ...(6) 
Let c/2 and d/2 be the heights of a line spectrum pair at 
25 mirror-image positions after filtering. 

The sum of sine waves at this time is expressed like 
equation (5) : 
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g(x) = d X cos(2 7rfgx) 

+ c X cos(27c(fs - fg)x) 
= 2d X cos(27c(fg/2 - fs)x) 
X cos (2 7c f sx/2) 

+ (c - d) X cos(27c(fs - fg)x) ...(7) 
The first term of equation (7) represents beat 

component; and the second term, a general sine wave 

component . 

According to equation (5), we have 



sin^ n ^ sin 7C — 

10 c / 2 = (a / 2) — -A ) ,d / 2 = (a / 2) ^ V ^ / (f g > fs) 



... (8) 

If in equation (7) the power of the second term which 
represents the normal sine wave component exceeds the power 
of the first term which represents the beat component, it 
15 may become difficult for the observer to recognize the beat 
component . 

If the ratio of the power of the first term to that 
of the second term in equation (7) is calculated on the basis 
of the above assumption, then 
20 (c - d)^/2d^ > 1 (condition under which the power of 

the second terminal exceeds that of the first term) 

This inequality can be rewritten into 

c/d > 2^^^ +1 

A substitution of this into equation (7) yields 

25 
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f s - fg 



fg 



> ^/^^TT, f g > 



1 



fs 



+ 1 



fg > 0.608 fs 



. . . (9) 



It is obvious from inequality (9) that if the 
frequency of a sine wave to be sampled is higher than 60.8% 
of a sampling frequency (aliasing occurs) , sampling can be 



grid stripe pattern information appears as a component 
equal to or less than 80% of the Nyquist frequency (fs/2) . 
This component is equivalent to a component equal to or less 
than 40% of the sampling frequency. If, therefore, the grid 
stripe pattern has a frequency equal to or less than 80% 
of the Nyquist frequency, a stable stripe pattern can be 
observed without any conspicuous beat upon sampling. 

The above consideration defines the upper limit 
spatial frequency of a stripe pattern (grid stripe pattern) 
appearing below the Nyquist frequency due to the grid. More 
specifically, the frequency of the grid stripe pattern is 
set to be equal to or less than 80% (equal to or less than 
40% of the sampling frequency) of the Nyquist frequency. 
In practice, however, this frequency has its own lower limit . 
Since an artifact originating from the stripe pattern 



performed with little observable beat noise. In this case. 
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itself, i.e., stripe pattern information, is not allowed 
to overlap an image component, the maximum frequency of the 
image component needs to be lower than the lower limit 
frequency of the grid stripe pattern. 
5 In general, the maximum frequency component of a 

signal representing an image cannot be accurately defined. 
Examples of evaluation criteria for images will be listed 
below: 

• A frequency that meets the approval of the observer 
10 when an acquired image is reproduced by a display apparatus 

or recording apparatus. 

• When a maximum frequency is assumed, the sampling 
pitch is determined by regarding 1.5 to 2 times the assumed 
frequency as a maximum frequency (see Nakamizo et al., " 

15 Counting/Measurement" (Baihukan) ) . 

The latter condition is widely used, in particular, 
and the sampling pitch is often set such that a spatial 
frequency required generally is equal to or less than 60% 
of the Nyquist frequency (equal to or less than 30% of the 

20 sampling frequency) . That is, the sampling pitch is set 
such that grid stripe pattern information appears at a 
frequency equal to or higher than the sampling frequency. 

Note that the lower limit of grid stripe frequencies 
is 25% the sampling frequency with reference to ^'2 times" 

25 of "1.5 to 2 times" described above and 1/3 the sampling 
frequency with reference to ^'1.5 times'' thereof. 
Therefore, the lower limit of grid stripe frequencies may 
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be set to at least 25% of the sampling frequency, preferably 
30% thereof/ more preferably 1/3 thereof. 

Assume that the sampling pitch is 0.1 mm (fs = 10 
cyc/mm) . In this case, the Nyquist frequency is 5 cyc/mm. 
5 According to the above general condition, the frequency of 
an image component which is generally used is equal to or 
less than 30% of the sampling frequency, i.e., 3 cyc/mmi. 
That is, the lower limit frequency of stripe pattern 
information is 3 cyc/mm, and the upper limit frequency is 
10 40% of the sampling frequency, i.e., 4 cyc/mm. If, 

therefore, the frequency of a grid stripe pattern is set 
within the range of 3 to 4 cyc/mm, no conspicuous beat 
appears in a grid stripe pattern and interferes with 
observation. 

15 The above condition is set to determine the frequency 

of a grid stripe pattern, and the frequency of the lead 
members of the grid body whose stripe pattern corresponds 
to this value changes after sampling. In this case, as 
shown in Fig, 9 a frequency range RG [cyc/mm] of the grid 

20 body is given by 

5(2n + 1) - 2 - 5(2n + 1) - 1 or 5(2n + 1) + 1 - (2n 
+ 1) + 2; n = 0, 1, 2,... (see Fig. 9) 
where n is an integer equal to or more than 0 . 

This mathematical expression represents a value 

25 calculated when the sampling pitch is 0 . 1 mm ( f s = 10 cyc/mm) . 
In general, if the sampling frequency fs [cyc/mm] is 
(sampling pitch 1/fs [mm] ) , the frequency range of the grid 
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body is given by 

fs(n +0.3) ~ fs(n + 0.4) 

or 

fs(n + 0.6) - fs(n + 0.8) [cyc/mm] 
5 The frequency of the grid body (the number of grid 

elements) is determined within the above range in 
consideration of scattered ray removing performance as the 
primary object of the grid, the resolution of a flat panel 
sensor, and the like. In general, in consideration of the 

10 resolution of the sensor, a grid stripe pattern frequency 
of 6 to 7 cyc/mm can be selected, at which the second-order 
harmonic is difficult to resolve and the scattered ray 
removing ratio is high. 

In this embodiment, the frequency of the grid to be 

15 used is selected by the above calculation such that the 
observer can be satisfied to some degree without removing 
any grid stripe pattern information owing to some 
experience or depending on an application purpose, thus 
overcoming the problem. 

20 As the spatial frequency of the grid stripe pattern 

is fixed, the grid stripe pattern information can be removed 
to some extent by filtering. In removing the grid stripe 
pattern information by setting a grid stripe pattern 
frequency in the above manner, even if the pattern 

25 information cannot be completely removed, reducing the 
intensity of the grid stripe pattern will minimize the 
influence on the observer. 
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Fig. 1 is a schematic view of the first embodiment - 
Fig. 1 shows a .system for imaging the human body lying on 
a table. Reference numeral 1 denotes an X-ray tube; 2, a 
human body as an object to be imaged; 11, a grid for removing 
5 scattered X-rays, which is a detachable grid for removing 
scattered X-rays from the object 2; 3, an X-ray sensor panel 
for converting an X-ray intensity distribution (X-ray 
transmission distribution) into a charge distribution, 
two-dimensionally sampling the distribution at desired 

10 intervals, and sequentially outputting the sampled data; 
5, an analog/digital converter; 4, a controller f or • 
controlling the X-ray radiation timing and image 
acquisition timing; and 6, a memory for temporarily storing 
an image. The X-ray sensor panel 3 varies in offset and 

15 gain for each pixel. To correct this variation, an offset 
value as an image acquired without any radiation of X-rays 
is stored in a memory 8, whereas data obtained by 
logarithmically converting a gain value acquired without 
the object 2 and grid 11 is stored in a memory 9. Reference 

20 numeral 7 denotes a conversion unit for logarithmic 

conversion, and more specifically, a lookup table. An 
acquired image of the human body is logarithmically 
converted after the offset value in the memory 8 is 
subtracted (removed) from the image. The difference 

25 between the resultant value and the gain value in the memory 
9 is calculated (division) to obtain an X-ray intensity 
distribution image having undergone correction of a 
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variation in gain value. This image is temporarily stored 
in a memory 10. Thereafter, the stored image is extracted 
and subjected to image storage, image processing, image 
display, and hardcopy operation, and the like to be used 
5 for diagnosis and the like, 

A block 12 is an image processing means (filtering 
means) for removing grid stripe pattern information by 
image processing (filtering) . The image processing means 
12 removes a grid stripe pattern component by spatial 

10 filtering using the image stored in the memory 10. A 

mechanism 13 is a selection means (switch) that is operated 
by the operator to choose between using the filtering means 
12 or not using it in accordance with an output from an 
external output from a grid removing operation effective 

15 setting means (mechanism) 16 or operation panel 20. The 
flow of a signal can be changed to skip the operation of 
the block 12 in accordance with the selection made through 
the selection means 13. 

In this case, the X-ray sensor panel 3 has a plurality 

20 of pixels distributed vertically and horizontally at a 
pitch of 0.1 mm in a two-dimensional space. With this 
structure, two-dimensional, discrete sampling is performed. 
As described in relation to the setting of the above grid 
stripe pattern frequency, the grid stripe pattern frequency 

25 (the number of grid elements) of the grid 11 is set to 6 
to 7 cyc/mm, and hence the observer can observe an image 
without a sense of incongruity even if the grid stripe 
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pattern information is not removed by filtering. 

Fig. 2A schematically shows the state of a 
one-dimensional amplitude spectrum in a direction 
perpendicular to the grid stripe pattern of the image stored 
5 in the memory 10. Referring to Fig. 2A, reference numeral 
32 denotes a spectrum of an image component; and 31^ a 
spectrum of a grid stripe pattern component, which exhibits 
a substantial spectrum form with noise being neglected. 
Even if the grid stripe pattern component 31 exists, the 

10 observer observes only this stable frequency component, and 
there is no component associated with beat noise. 
Therefore, as the observer becomes accustomed to the stripe 
pattern or recognizes it, the existence of the stripe 
pattern does not relatively interfere with the observation. 

15 If, however, the operator or observer wants to remove 

the grid stripe pattern information from this image owing 
to subsequent image processing or the image reproducing 
mechanism, he/she selects grid removing operation through 
the operation panel 20. 

20 Fig. 2B schematically shows the state of a spectrum 

upon execution of grid removing filtering. Reference 
numeral 33 denotes an example of a filter characteristic; 
and 34, an image spectrum after filtering. To stabilize 
spatial characteristics, the filter cannot have a steep 

25 characteristic. If, therefore, the grid stripe pattern 
component 31 is removed, part of an image component 
inevitably deteriorates in response characteristic. In 
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consideration of this, the operator chooses between removal 
and nonremoval of the grid. In this embodiment, when the 
acquired image exhibits the gain value stored in the memory 
9, no problem arises if the image is acquired without 
5 removing the grid 11 because the frequency of the grid is 
constant • 

Fig. 3 is a flow chart for the implementation of this 
embodiment by means of software. Referring to Fig. 3, a 
process block (step) is divided into operations in blocks 

10 CI to Cll. The operation in block CI is executed to acquire 
an image of a gain value. An image A is obtained by 
radiating X-rays without any object. In block C2, this 
image is logarithmically converted into an image B. In 
block C3, an offset value is obtained; an image C is obtained 

15 without radiating any X-rays. In block C4, an image of the 
object is actually acquired; the grid is installed, and an 
image D is obtained by irradiating the object with X-rays. 
In block C5, the image C is subtracted from the image D to 
obtain an image E having undergone offset correction. In 

20 block C6, the image E is logarithmically converted into an 
image F. In block CI, the image B is subtracted from the 
image F to obtain an object image G having undergone gain 
correction. In block C8, the flow branches depending on 
whether a stripe pattern operating from grid will be removed 

25 in accordance with operation (instruction) of the operation 
panel 20 by the operator. In block C9 , since the 
instruction to remove the grid is received, filtering, i.e.. 
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grid removing operation, is performed for the image G to 
obtain an object image H from which grid stripe pattern 
information is removed. In block CIO, the object image H 
is output. If no grid removing instruction is received, 
5 the object image G is output without any processing in block 
Cll. 

In this embodiment, as the memory means 10 in Fig. 1, 
a nonvolatile storage medium such as a magnetic disk may 
be used to always store image data including grid stripe 

10 pattern information to allow the operator to select an image 
without any grid stripe pattern information or an image with 
grid stripe pattern information or output them at once. 
Furthermore, this embodiment can be practiced with an 
arrangement including only an image storage system without 

15 any image acquisition system. 
(Second Embodiment) 

Fig. 4 is a block diagram showing the second 
embodiment, in which the frequency of a grid stripe pattern, 
the sampling pitch of an X-ray sensor panel, and the like 

20 are set in the same manner as in Fig. 1. In the structure 
shown in Fig. 4, an application purpose table 14 is prepared. 
After an image is acquired, the operator selects an 
application purpose through an application purpose setting 
means 17 with respect to the image stored in a memory 10 

25 or magnetic disk. As a consequence, whether to remove the 
grid or not is automatically selected. If a switch 13 
selects the A side, the image in the memory 10 is filtered 
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by a filtering unit 12 and output. If the switch 13 selects 
the B side, the image in the memory 10 is output without 
being filtered. 

As described in relation to the setting of the above 
5 grid stripe pattern frequency, the grid stripe pattern 
frequency (the number of grid elements) of a grid 11 is set 
to 6 to 7 cyc/mm, and hence the observer can observe an image 
without a sense of incongruity even if the grid stripe 
pattern information is not removed by filtering. 

10 With this table 14 , for example, in emphasizing a high 

spatial frequency as in spatial frequency emphasis 
processing^ since a grid .image becomes a hindrance, it is 
removed (the switch 13 is set on the A side) . In displaying 
an image or performing hardcopy operation on a larger scale, 

15 i.e., 100% or more, the switch 13 is set on the B side to 
inhibit removal of a grid image so as to minimize an image 
blur. In displaying an image or performing hardcopy 
operation upon reduction, the switch 13 is set on the A side 
to remove a grid image. If an image is to be stored in 

20 another storage means, since removal processing can be 
performed for the stored image, the switch 13 is set on the 
B side to inhibit removal of a grid image, thus increasing 
the information amount. 
(Third Embodiment) 

25 Fig. 5 is a block diagram showing the third 

embodiment, in which the frequency of a grid stripe pattern, 
the sampling pitch of an X-ray sensor panel, and the like 
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are set in the same manner as in Fig. 1. In the structure 
shown in Fig. 5, an imaging position table 18 is prepared. 
When the operator selects an application purpose through 
an imaging position setting means 19 with respect to the 
5 image stored in a memory 10 or magnetic dick upon imaging 
operation, whether to perform grid removing operation or 
not is automatically selected. 

As described in relation to the setting of the above 
grid stripe pattern frequency, the grid stripe pattern 

10 frequency (the number of grid elements) of a grid 11 is set 
to 6 to 7 cyc/mm, and hence the observer can observe an image 
without a sense of incongruity even if the grid stripe 
pattern information is not removed by filtering. 

With this table 18, when the observer is to observe 

15 an image of a bone portion such as the pelvis or joint, which 
requires a high spatial frequency for the image, a switch 
13 is set on the B side to allow the observer to observe 
the image without any blur without removing a grid image. 
When the observer is to observe a chest portion (front chest 

20 portion) , abdomen, or the like for which a high spatial 
frequency is not required, and a halftone image needs to 
be easily observed, the switch 13 is set on the A side to 
allow the observer to observe the image without any grid 
image . 

25 (Fourth Embodiment) 

Although the grid frequency is so set that a grid image 
does not relatively interfere with observation, if the 
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contrast of the grid image is strong, it still interferes 
with the observer. The contrast of a grid image varies 
depending on the conditions (e.g., energy) of X-rays to be 
used. In some case, an object is imaged without any grid. 
5 To determine this, the spectrum of a given portion 

of an acquired image in a direction perpendicular to the 
grid is calculated, and whether to perform grid removing 
operation or not is selected depending on the peak value 
of the grid component (directly corresponding to the 

10 contrast value if it is a logarithmic image) . 

Fig. 6 is a block diagram showing the fourth 
embodiment, in which the frequency of a grid stripe pattern, 
the sampling pitch of an X-ray sensor panel, and the like 
are set in the same manner as in Fig. 1. As described in 

15 relation to the setting of the above grid stripe pattern 
frequency, the grid stripe pattern frequency (the number 
of grid elements) of a grid 11 is set to 6 to 7 cyc/mm, and 
hence the observer can observe an image without a sense of 
incongruity even if the grid stripe pattern information is 

20 not removed by filtering. 

A block 15 in Fig. 6 includes software; a flow chart 
is shown in the block. In block C21 in the flow chart, an 
arbitrary line is acquired from an image from a memory 10. 
In block C22, one-dimensional Fourier transformation is 

25 performed for this image to calculate an amplitude spectrum. 
In block C23, an amplitude value (spectrum value) Vp 
corresponding to the grid frequency is measured from the 
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amplitude spectrum value. In block C24, the amplitude 
value Vp is compared with a threshold TH set in advance. 
If Vp is larger than TH, the corresponding grid image must 
be removed. To remove the grid image, therefore, a switch 
5 13 is set on the A side in block C25. If Vp is not larger 
than TH, the switch 13 is set on the B side in block C26. 

Note that the switch 13 may be selectively operated 
in accordance with the magnitude of contrast of a grid 
stripe pattern existing at the grid frequency. If the 
10 contrast is higher than a predetermined threshold, the 
switch 13 is set on the A side to remove a grid image. If 
the contrast is lower than the predetermined threshold, the 
switch 13 is set on the B side to inhibit the removal of 
a grid image. 

15 As described above, according to this embodiment, in 

a system for two-dimensionally sampling an X-ray image and 
forming a digital image, an image that can be observed by 
the observer with little sense of incongruity without 
removing a grid stripe pattern from the image can be formed 

20 by setting the frequency of the grid lead members to be equal 
to or lower than 40% of the sampling frequency at which the 
spatial frequency of a stripe image originating from the 
grid for removing scattered rays from the image and equal 
to or higher than the frequency (60% of the Nyquist 

25 frequency in general) at which a grid image does- not easily 
interfere with observation of the image by the observer and 
does not overlap an image component . 
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In addition, the filtering means 12 can choose 
between removal or nonremoval of a grid stripe pattern in 
accordance with selection through the switch 13. Since a 
grid stripe pattern is automatically removed in accordance 
5 with operation by the operator (observer) , the application 
purpose of an image, imaging position, and the amplitude 
(intensity) of a grid stripe pattern, an appropriate image 
can be obtained by removing a grid stripe pattern only when 
required. 

10 Each embodiment described above is merely an example 

in practicing the present invention. Note that the 
technical scope of the present invention should not be 
interpreted in a limited manner. That is, the present 
invention can be practiced in various forms without 

15 departing from the spirit and scope of the present 
invention. 

As has been described above, according to each 
embodiment described above, when an X-ray image is 
two-dimensionally sampled to acquire a digital image, an 

20 image that can be observed by the observer with little sense 
of incongruity without removing a grid stripe pattern from 
the image can be formed by setting the spatial frequency 
of a stripe pattern originating from the grid for removing 
scattered rays from the image to a predetermined value. 

25 In addition, since whether to remove a stripe pattern 

originating from the scattered ray removing grid from an 
image can be selected, an appropriate image can be obtained 
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by removing a stripe pattern originating from the scattered 
ray removing grid. 
(Fifth Embodiment) 

In the above case, since the resolving power of the 
5 X-ray image receiving unit (X-ray image sensor) is not very 
high, the above apparatus is based on the premise that the 
grid pattern exhibits an almost single sine wave. This will 
be described below • 

Recently, with technical innovation of X-ray image 

10 sensors, a scheme of directly converting an X-ray intensity 
distribution into a charge distribution by using, for 
example, a method of acquiring free electrons, generated 
by X-rays, using a strong electric field, has been studied 
and put into practice, in place of an indirect conversion 

15 scheme of converting an X-ray intensity distribution into 
a fluorescence distribution first and then 
photoelectrically converting the fluorescence 
distribution. 

In the case of the direct conversion scheme, an 

20 aperture used to acquire X-rays as charges (electrons) is 
the main factor that determines resolving power, and hence 
a high-resolving power X-ray image sensor can be obtained. 
When such a high-resolving power X-ray image sensor is used, 
an original grid pattern is resolved more finely (harmonic 

25 components are also resolved) regardless whether the sensor 
is not based on the direct conversion scheme. It is 
therefore expected that a grid pattern in an acquired image 
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may not exhibit a single sine wave. 

Such a situation will be described with reference to 
Fig. 15. Fig. 15 is a graph on a one-dimensional spatial 
frequency axis, with the abscissa representing the spatial 
frequency. Referring to Fig. 15, a Nyquist frequency Fn 
in the middle is a Nyquist frequency set when the pixels 
of the sensor are regarded as pixels corresponding to 
spatial sampling, i.e., 1/2 a sampling frequency (the 
reciprocal of the sampling pitch) Fs. 

The peak denoted by reference numeral 201 in Fig. 15 
is the spatial frequency of the lead foil members of the grid, 
which is set to Fg = 1.25Fn, for the sake of convenience. 
According to the sampling theorem, this frequency is 
expressed by a frequency equal to or lower than the Nyquist 
frequency, and the frequency Fml can be calculated as 
follows : 

Fml = 2-Fn - Fg = 0.75Fn ...(10) 
This peak is denoted by reference numeral 202 in 
Fig. 15 . If the resolving power of the sensor in use is high, 
a second-order harmonic 2Fg of a grid frequency Fg is sampled 
at the same time. A frequency Fm2 of this harmonic is 
calculated as follows and is denoted by reference numeral 
203 in Fig. 15: 

Fm2 = 2-Fg - 2-Fn = 0.5Fn ...(11) 
A third-order harmonic Fm3 is denoted by reference 
numeral 204 in Fig. 15 and given by 

Fm3 = 4-Fn - 3-Fg = 0.25Fn ...(12) 
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In general, a frequency Fk on an image with respect 
to k-th order harmonic is expressed by 

Fk = |2- j-Fn - k-Fg| . . . (13) 

(An integer j including 0 is selected to satisfy 0 ^ Fk 
5 ^ Fn.) 

As described above, when a high-resolving power sensor 
is used, even line spectra originating from such harmonics 
appear as stripe pattern information on an image. According 
to the conventional technique, in order to remove such stripe 

10 pattern information originating from a plurality of spatial 
frequencies from an image, filtering corresponding to each 
frequency must be done. This inevitably affects image 
information. 

The following embodiment has been made in 

15 consideration of the above problems and explains a radiation 
imaging method and apparatus which facilitate removing or 
reducing components originating from a grid in an acquired 
image or can acquire an image that is easy to observe even 
in the presence of such components, and a design method for 

20 the apparatus. 

In addition, the following embodiment explains a 
radiation imaging method and apparatus which facilitate 
removing or reducing harmonic components originating from 
a grid in an acquired image or can acquire an image that is 

25 easy to observe even in the presence of such components, and 
a design method for the apparatus. 

The following embodiment explains a radiation imaging 
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method and apparatus which can integrate a fundamental 
frequency component and second-order harmonic component 
originating from a grid in an acquired image into a 
substantially single spectrum, and a design method for the 
5 apparatus. 

Furthermore, the following embodiment explains a 
radiation imaging method and apparatus which can integrate 
components originating from a grid in an acquired image into 
a substantially single spectrum, and a design method for the 
10 apparatus. 

Other objects of the following embodiment will be 
obvious from the following description of the specification. 

According to the following embodiment, image 
components originating from a grid in an acquired 
15 (visualized) image are integrated (converged) in a 
substantially single spectrum. 

The fifth embodiment will be described below with 
reference to the accompanying drawings. 

The fifth embodiment exemplifies an arrangement for 
20 converging grid stripe spectra up to at least the 

second-order harmonic into a substantially single spatial 
frequency spectrum by spatial sampling. This arrangement 
makes a grid stripe pattern (an image component originating 
from the grid) become an approximate sine wave form and 
25 prevent the occurrence of unnecessary low-frequency 

components. In addition, a grid pattern can be easily or 
properly removed, as needed, by grid pattern removal 
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processing with subsequent filtering operation or the like. 

Conditions under which the fundamental frequency (k 

=1) of a grid pattern component coincides with the 

frequency (k = 2) of the second-order harmonic component 

in a sampled image are calculating by referring to equation 

(13) . Since this is an absolute value computation, the 

following two conditions will be considered. 

2jiF, - = 2j,F, - 2Fg ... (14 - 1) 

2jiF„ - F, = 2F, - 2J,F, . . . (14 - 2) 

where ji and j2 are positive integers (natural numbers) 
including 0. 

By solving equation (14-1) for Fg, the following 
equation is obtained: 

Fg = 2Fn(j2 - ji) ... (15) 

This coincides with the condition described in 
Japanese Patent Laid-Open No. 9-75332. If Fg slightly 
deviates from this condition due to a manufacturing error 
in the grid or the like, a stripe pattern with a very low 
frequency is produced. This greatly damages image 
information and makes it difficult to remove a grid stripe 
pattern by filtering or the like. 

By solving equation (14-2) for Fg, the following 
equation is obtained: 

Fg = + j2)Fn = ^(^1 . • . (16) 

In this case, Fs = 2-Fn (sampling frequency) . 

In this case, (ji -f 32) is a positive integer (natural 
number) and can take numerical values such as (2/3) Fn, 
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(4/3) Fn, and (6/3) Fn. 

If, however, (ji + 32) is a multiple of three, the 
resultant condition is the same as that represented by 
equation (15) . Since this condition is not appropriate as 
5 described above, this case is excluded. If the grid 

frequency Fg is set to satisfy equation (16) when (ji + J2) 
is a natural number excluding a multiple of three, there 
is no possibility that a low-frequency stripe pattern like 
that described above will be produced, even with a slight 

10 manufacturing error. 

In addition, if the grid frequency is set in 
accordance with this condition, the fundamental frequency 
of the grid stripe pattern and the frequency of the 
second-order harmonic become almost equal to each other in 

15 the image (become a substantially single spectrum) , and 
only a substantially single stripe pattern with a sine wave 
is advantageously generated in the image. The present 
inventor has found that when the second-order harmonic 
component of a grid stripe pattern in an acquired image 

20 cannot be neglected, a very good effect can be obtained by 
setting the grid frequency so as to satisfy equation (16) . 

In the case of general medical images, the sampling 
pitch is about 100 /zm to 200 //m, and hence the Nyquist 
frequency Fn falls within the range of about 2.5 cyc/mm to 

25 5 cyc/mm. The frequency Fg of a physically appropriate grid 
as a grid aimed at removing scattered X-rays falls within 
the range of about 3 cyc/mm to 10 cyc/mm. For this reason, 
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a 



( ji + j2) =2 exhibits a promising condition that may coincide 
with the condition represented by equation (16) . The 
condition represented by (ji + j2) = 3 coincides with the 
condition represented by equation (15), and hence is 
5 inappropriate. Therefore the condition represented by 
equation (17) is appropriate. 




• (17) 



In this embodiment, by making the grid frequency 
almost coincide with a frequency that satisfies the 

10 condition represented by equation (16), the spectra of an 
essentially unnecessary grid stripe pattern in a sampled 
image can be intefrated to near Fs/3 (= 2Fn/3) . This makes 
it easy to remove a grid stripe pattern from an image by 
filtering or the like. 

15 According to the foregoing embodiments (Japanese 

Patent Application No. 2000-028161 filed by the present 
applicant), the grid frequency Fg is set such that the 
frequency of an image component originating from a grid in 
acquired image data becomes equal to a frequency selected 

20 from the range from 50% the Nyquist frequency Fn to 80% 
thereof. This condition makes moire fringes (beat-like 
variation components) originating from a grid less 
noticeable even if imaging is done with the fixed grid, and 
includes the condition Fg = (2/3) jFn (j is a natural number 

25 excluding a multiple of three) according to this 
embodiments . 



The solution (condition) in this embodiment which is 
indicated by equation (16) is a special solution obtained 
by also taking harmonic components (second-order harmonic 
component) of a grid stripe pattern into consideration, and 
5 is also made with respect to the technical problem newly 
found by the present inventor. This solution has the unique 
effect described above. Therefore, the invention 
described in this embodiment constitutes valid selection 
inventions and numerically limited inventions. 

10 Fig. 11 is a block diagram schematically showing an 

X-ray image acquisition apparatus according to the fifth 
embodiment of the present invention • A case wherein the 
present invention is applied to an X-ray image acquisition 
apparatus for imaging an object to be imaged by using X-rays 

15 will be described below. However, the present invention 
can also applied to an image acquisition apparatus using 
radiation other than X-rays. Referring to Fig. 11, 
reference numeral 101 denotes an X-ray generating unit 
which has an X-ray tube, high voltage generator, and 

20 controller and emits X-rays in the direction indicated by 
the arrow; 102, an object to be imaged, typically the human 
body; 103, a support base such as a bed which supports a 
lying object such as the human body; and 104, a scattered 
ray removing grid having lead foil members whose intervals 

25 are optimally set. This setting will be described later. 

Reference 105 denotes an X-ray image sensor which 
converts the intensity distribution of X-rays (X-ray image) 
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transmitted through the object into an electric signal and 
is formed by using a large-size solid-stage image sensing 
element having a detection area (image-receiving area) 
constituted by a plurality of pixels arranged 
two-dimensionally in the form of a matrix. The X-ray image 
sensor 105 spatially samples an X-ray image through the 
detection area. 

In this embodiment, as this X-ray image sensor, a 
sensor exhibiting little decrease in resolving power in an 
energy conversion process (e.g., a sensor based on the 
direct conversion scheme) is used. In addition, the 
sampling pitch (pixel pitch) of this X-ray image sensor is 
set to 0.16 mm. The X-ray image sensor will be referred 
to as a flat panel sensor or simply as a sensor hereinafter. 
The flat panel sensor is controlled by a controller (not 
shown) to sequentially scan charges proportional to the 
amounts of X-rays present for the respective pixels and 
convert them into predetermined electrical quantities 
(e.g., voltages or currents), thereby outputting X-ray 
image information as electrical, quantities. 

Reference numeral 106 denotes an A/D converter for 
converting the analog electrical quantities output from the 
flat panel sensor 105 into digital values; 107, a memory 
(storage unit) for temporarily storing a digital value from 
an A/D converter as image information; 108, a switch for 
reading out the information stored in the memory 107 and 
switching destinations for the information; 109, a memory 
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(storage unit) for storing the image signal output from the 
X-ray image sensor 105 without irradiation of X-rays as an 
offset fixed pattern image; and 110, a memory (storage unit) 
for storing an image of the object which is obtained by 
5 actually irradiating the object with X-rays. 

According to a specific imaging method, an X-ray dose 
measuring device (not shown) called a phototimer which 
monitors the amount of X-rays transmitted through an object 
to be imaged is used to control X-ray emission by the X-ray 

10 generating unit, and X-ray emission by the X-ray generating 
unit is stopped at the instant when the cumulative amount 
of X-rays transmitted through the object becomes a 
predetermined value. At the same time when X-ray emission 
is stopped, the controller for this X-ray image acquisition 

15 apparatus causes the flat panel sensor to start scanning, 
and temporarily stores image information of the object 102 
in the memory 107. Thereafter, the controller sets the 
switch 108 to the A side to make the memory 110 store the 
image information. Immediately after this operation, the 

20 controller causes the flat panel sensor to store charges 
without emitting X-rays for the same time as the X-ray 
emission time determined by the above phototimer, and then 
makes the flat panel sensor scan. The controller then 
causes the memory 107 to temporarily store the resultant 

25 image information as an offset fixed pattern image. After 
this operation, the controller sets the switch 108 to the 
B side to make the memory 109 store the offset fixed pattern 



image from the memory 107. Reference numeral 111 denotes 
a subtracter which sequentially subtracts, essentially 
from the values of object image data in the memory 110, the 
values of the offset fixed pattern image data in the memory 
109 at the corresponding positions; and 112, a memory 
(storage unit) which stores the difference (image data 
after offset fixed pattern correction) from the subtracter 
111. 

Reference numeral 113 denotes an LUT (Look Up Table) 
for converting image data into a logarithmic value; 114, 
a switch for switching destinations for data from the LUT 
113; and 115, a memory (storage unit) for storing data from 
the switch 114 whose output destination is set to the C side. 
In general, image data from the memory 112 which is obtained 
through an object to be imaged is logarithmically converted 
by the LUT 113. The image data is then stored in the memory 
115 through the switch 114 set to the C side. 

Reference numeral 116 denotes a memory (storage unit) 
for storing image data when this X-ray image acquisition 
apparatus performs operation called calibration imaging. 
In this calibration imaging, imaging is performed in the 
same manner as described above, and the switch 114 is set 
to the D side. The resultant image data is then stored in 
the memory 116. This operation differs from the above 
object imaging operation in that imaging is performed by 
emitting X-rays without the mediacy of the object 102 . With 
this operation, image data (simply referred to as gain 
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variations or gain variation data) is stored in the memory 
116, which constitutes gain variations (also called 
sensitivity variations) of a plurality of pixels 
constituting the flat panel sensor and the intensity 
5 distribution (shading) of X-rays emitted from the X-ray 
generating unit. In general, this calibration imaging is 
done about once every day at the start of work. 

Reference numeral 117 denotes a gain corrector 
(subtracter) having the function of correcting image data 

10 on the basis of gain variations by subtracting the gain 
variation data stored in the memory 116 from the image data 
stored in the memory 115 for each corresponding pixel (this 
operation substantially equivalent to division because 
each data having undergone logarithmic conversion) ; and 118, 

15 a memory (storage unit) which stores object image data 
having undergone correction based on gain variations which 
is performed by the gain corrector 117. 

Grid frequency (lead foil member pitch) selection 
will be described below. Fig. 12 schematically shows the 

20 flat panel sensor 105 according to this embodiment. Each 
portion 300 indicated by a rectangle represents a portion 
for receiving X-ray energy (charge) upon conversion, i.e., 
a pixel. As described above, the pixel pitch of this sensor 
is set to 0.16 mm. It is, however, difficult for reasons 

25 of manufacturing techniques to match the size of one side 
of a pixel, i.e., aperture, with the pixel pitch, and hence 
the size is set to 0.14 mmD . As described above, since 
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the flat panel sensor 105 is an X-ray image sensor of the 
type that exhibits little decrease in rresolving power in 
an energy conversion process, the MTF (Modulation Transfer 
Function) of the sensor is mostly determined by the shape 
of this aperture. 

Referring to Fig. 13, reference numeral 301 denotes 
the MTF of the the flat panel sensor 105. In this graph, 
the ordinate represents the MTF; and the abscissa, the 
spatial frequency. This MTF corresponds to the data 
obtained by Fourier transformation of an aperture of 0.14 
mm. Reference numeral 305 denotes the sampling frequency 
of the flat panel sensor 105, which is 1/0.16= 6.25 cyc/mm; 
and 306, a Nyquist frequency which is 6.25/2 = 3 . 125 cyc/mm. 
In this case, the pitch of the shades of the lead foil members 
of the grid on the flat panel sensor is set to 0.24 mm. As 
a consequence, the fundamental spatial frequency of the 
grid stripe pattern becomes 1/0.24 =4.17 cyc/mm. 
Referring to Fig. 13, this frequency is denoted by 
reference numeral 303. The grid stripe pattern is the shade 
pattern of the lead foil members. Strictly speaking, 
therefore, the shade pattern cannot be a sine wave but has 
a plurality of harmonics, and the frequency spectra of the 
shade pattern are a set of a plurality of line spectrum. 
Reference numeral 304 denotes the second-order harmonic of 
the grid stripe pattern, which has a frequency of ( 1/0 . 24 ) *2 
= 8.33 cyc/mm. 

As also described with reference to equation (17), 
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if the grid strip pattern is visualized below the Nyquist 
frequency by sampling, the grid stripe patterns 303 and 304 
are integrated into a single spectrum (line spectrum) with 
a frequency of 2 . 08 cyc/mm. In this case and general cases, 
5 since third- or higher-order harmonics are hardly resolved 
even by this MTF defined by only the aperture, this spectrum 
is the only spectrum (line spectrum) of the grid stripe 
pattern which appears in an image. 

Referring back to Fig. 11, reference numeral 119 

10 denotes an image processing unit which performs 
one-dimensional spatial filtering in a direction 
perpendicular to the grid stripe pattern. The image 
processing unit 119 selectively removes the grid stripe 
pattern 302 in Fig. 13. More specifically, this removal 

15 can be realized as follows. A grid stripe pattern is 

estimated and formed such that image components including 
the grid stripe pattern extracted from image dada by 
filtering are processed on the basis of the essential 
characteristics of the grid which exhibit stable 

20 periodicity. The formed grid stripe pattern is then 

subtracted from the original image data (this operation is 
substantially division because the data have undergone 
logarithmic conversion) (this method is disclosed in 
Japanese Patent Application No. 2001-134208 filed by the 

25 present applicant) . Reference numeral 120 denotes a 
memory (storage unit) which stores image data having 
undergone processing by the image processing unit 119. 
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Note that integration of the fundamental wave and 
second-order harmonic wave of the grid stripe pattern into 
substantially one spatial frequency as in this embodiment 
is suited to easily or properly removing the grid stripe 
5 pattern regardless of the grid stripe pattern removing 
method to be used. 

Fig. 14 shows the positional relationship between an 
X-ray source (X-ray focal point) , a grid, and a flat panel 
sensor. Reference numeral 4 01 denotes an X-ray source; 4 02, 

10 a grid; and 403, a flat panel sensor. Depending on the 
target region of an object to be imaged, the amount of 
scattered X-rays is small, and imaging is performed without 
using the grid. For this reason, an X-ray image acquisition 
apparatus is generally designed to selectively use a grid. 

15 For example, a detachable grid is mechanically mounted in 
front of the sensor. Consequently, a small gap D is spaced 
between the grid and the flat panel sensor, as shown in 
Fig. 14. In general, an X-ray source is a nearly point ray 
source and placed at a distance L from the grid. The 

20 plurality of lead foil members of the grid are not simply 
arranged parallel to each other but are arranged at angles 
that make them converge to the X-ray source at this distance . 
In this case, since L and D are about 180 cm and 10 mm, 
respectively, a difference owing to a magnifying effect 

25 appears between the actual pitch of the lead foil members 
of the grid and the pitch of the shades of the lead foil 
members on the flat panel sensor. Letting gl be the actual 
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grid pitchy and g2 be the shade grid pitch on the sensor, 
the pitch g2 can be calculated by 

g2 = gl- (1 + D/L) ... (18) 

If the grid frequency (l/g2) on the sensor is set to 
5 4.17 cyc/mm as in the above case, the actual grid frequency 
(1/gl) is preferably set to 4.17- (1 + D/L) ^ 4.2 cyc/mm. 

A grid stripe pattern can be removed from the image 
data sampled by the flat panel sensor by one filtering 
operation, even if the spatial frequencies of the 
10 fundamental wave and second-order harmonic of the grid 
stripe pattern do not perfectly coincide with each other, 
as long as their difference falls within the range of errors 
of several %. Therefore, slight errors between them are 
allowed. This error ratio r is given by the following 
15 equation with reference to a target spatial frequency Fs/3: 

r = ELZjE3i ...,19, 

- Fs 
3 

In general, this error ratio preferably falls within 

5%. 

Equation (19) is generally expressed as 
20 ^ ^ t£s-_F^^Jj£s-JF2l ^^^^ 

- Fs 
3 

(Note that ji and j2 are so selected as to satisfy |jiFs - Fg| 
< Fs/2 and |j2Fs - 2Fg| < Fs/2.) 

Note that image processing unit 119 need not always 
remove a grid stripe pattern by filtering or the like unless 
25 the operator who observes the image requires so. According 
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to this embodiment, even if a grid stripe pattern is not 
removed by filtering or the like, since the grid stripe 
pattern in the image becomes a substantially single sine 
wave, it does not easily interfere with observation. 

This embodiment has exemplified the case wherein the 
ideal flat panel sensor exhibiting little decrease. in 
resolving power in an energy conversion process is used. 
However, the other arrangement of this embodiment can be 
effectively applied to a case wherein a flat panel sensor 
exhibiting a substantial decrease in resolving power in an 
energy conversion process is used. 

In addition^ all the grid frequencies Fg set when j 
= (ji + 32) in equation (16) takes natural numbers other 
than multiples of three can be used. 
(Sixth Embodiment) 

One of the characteristic features of a flat panel 
sensor is that pieces of information of a plurality of 
pixels adjacent to each other on the sensor can be added 
as analog values by accumulating the charges in the 
plurality of pixels on the sensor. Another characteristic 
feature is that the driving modes of the sensor can be 
switched, and the sensor can be designed to switch between 
the addition mode of performing such analog addition and 
the normal mode. 

When, for example, the dynamic function (movement) 
of the heart is to be observed, the flat panel sensor is 
driven at a high frame rate to acquire so-called moving 
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image data. In this case, for example, 2x2 adjacent pixels 
are handled together as one pixel, and image information 
is output from the flat panel sensor. If image data with 
higher resolution is required as a still image, the mode 
5 of the sensor is switched to the normal mode to increase 
the resolution. 

If the data of 2 x 2 pixels are added, the 
substantially sampling pitch is doubled. As a consequence, 
the Nyquist frequency Fn' becomes (1/2) Fn. Obviously, in 

10 the addition mode, with substitution of Fn = 2Fn' into 
equation (16), the grid frequency Fg satisfies the same 
condition as that in the normal mode. 

In the X-ray image acquisition apparatus according 
to this embodiment as well, if the grid frequency Fg is so 

15 set as to satisfy equation (16) in the normal mode (note 
that j = (ji + j2) is set to a natural number other than 
a multiple of three) , and the sampling frequency decrease 
ratio (which represents the number of pixels to be added 
in a direction perpendicular to the grid foil members, e.g., 

20 2 in the above case of addition of 2 x 2 pixels) due to pixel 
addition in the addition mode is set to a natural number 
other than a multiple of three, the spatial frequencies of 
the fundamental wave and second-order harmonic of the grid 
stripe pattern substantially coincide with each other in 

25 the image data sampled by the flat panel sensor regardless 
of the normal mode and addition mode. This makes it 
possible to obtain the same effect as that in the fifth 
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embodiment . 

As described above, by setting the grid frequency Fg 
so as to satisfy equation (16) (note that j = (ji + J2) is 
set to a natural number other than a multiple of three) , 
the frequencies of a grid stripe pattern in a visualized 
image can be converted (integrated) into a substantially 
single spectrum. This makes it possible to easily remove 
the grid stripe pattern from the visualized image. 

Note that the fifth and sixth embodiments are only 
examples in practicing the present invention. It should 
be noted that the technical range of the present invention 
is not limited by these embodiments. The present invention 
can be practiced in various forms without departing from 
the technical idea and main characteristic features of the 
present invention . 

As has been described above, according to the fifth 
and sixth embodiments, there are provided a radiation 
imaging method and apparatus which facilitate removing or 
reducing components originating from a grid in an acquired 
image or can acquire an image that is easy to observe even 
in the presence of such components, and a design method for 
the apparatus. 

In addition, according to the fifth and sixth 
embodiments, there are provided a radiation imaging method 
and apparatus which facilitate removing or reducing harmonic 
components originating from a grid in an acquired image or 
can acquire an image that is easy to observe even in the 
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presence of such components, and a design method for the 
apparatus. 

According to the fifth and sixth embodiments, there 
are provided a radiation imaging method and apparatus which 
5 can integrate a fundamental frequency component and 

second-order harmonic component originating from a grid in 
an acquired image into a substantially single spectrum, and 
a design method for the apparatus. 

Furthermore, according to the fifth and sixth 
10 embodiments, there are provided a radiation imaging method 
and apparatus which can integrate components originating 
from a grid in an acquired image into a substantially single 
spectrum, and a design method for the apparatus. 

The present invention is not limited to the above 
15 embodiments and various changes and modifications can be 
made within the spirit and scope of the present invention. 
Therefore, to apprise the public of the scope of the present 
invention the following claims are made. 
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